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We have searched for direct pair production of scalar top and scalar bottom quarks in 88 pb~^ 
of pp collisions at y^ =1.8 TeV with the CDF detector. We looked for events with a pair of heavy 



flavor jets and missing energy, consistent with scalar top quark decays to a charm quark and a 
neutralino, or scalar bottom quark decays to a bottom quark and a neutralino. The numbers of 
events that pass our selection for each process show no deviation from Standard Model expectations. 
We compare our results to next-to-leading order calculations for the scalar quark production cross 
sections to exclude regions in scalar quark-neutralino parameter space. 

PACS numbers: 14.80.Ly, 13.85.Rm 

Supersymmetry (SUSY) R] assigns to every fermionic Standard Model (SM) particle a bosonic superpartner and 
to every bosonic SM particle a fermionic superpartner. Therefore, the SM quark helicity states ql and qr acquire 
scalar partners qL and qR. SUSY models usually predict that the masses of the first two generations of scalar quarks 
are approximately degenerate. The scalar top quark (i) mass, however, may be lower than that of the other scalar 
quarks due to a substantial Yukawa coupling resulting from the large top quark mass. In addition, mixing between 
II and in can cause a large splitting between the mass eigenstates <i and ^2 0] • If the Yukawa coupling strength or 
the mixing is strong enough (or both) then the mass eigenstate fi can be lighter than the top quark. We note that 
many baryogenesis models require a light stop quark H . 

The bottom quark mass is much smaller than the top quark mass, therefore the effect of the Yukawa coupling on 
the scalar bottom quark (&) mass is small. However, in some regions of SUSY parameter space a large mixing between 
6l and ha can still occur, leading to a significant splitting between mass eigenstates and a low mass value for the 
lighter mass eigenstate (61) Q. 

At the Tevatron, third generation scalar quarks are expected to be produced in pairs via gg fusion and qq annihi- 
lation. In this Letter, we describe two analyses looking for processes in a Minimal Supersymmetric Standard Model 
framework: (i) a scalar top analysis, searching for the process pp — > titi -^ (cXi) (cXi)? and (ii) a scalar bottom 
analysis, searching for the process pp —> bibi —> (&Xi) (^Xi)- We do not assume a gaugino unification hypothesis. 
We assume the lightest neutralino Xi is the lightest supersymmetric particle and stable. This leads to experimental 
signatures with appreciable missing transverse energy. The decay ii -^ ex?, as in process (i), dominates via a one-loop 
diagram in the absence of flavor-changing neutral currents if nij: < irii, + 7ti-± and m^ < mw + fnb + "^v" H- For 

process (ii) we assume rrif^ > rrib + m^o and ttIj, < rrib + m^a y. Here, Xi and X2 ^^^ the lightest chargino and 
next-to- lightest neutralino. Therefore, the signature of both processes is a pair of acolinear heavy flavor jets, ^t, and 
no leptons in the final state. 

We have searched data corresponding to a total integrated luminosity of 88.0±3.6 pb~^ collected using the CDF 
detector during the 1994-95 Tevatron run. CDF is a general purpose detector and is described in detail elsewhere [||. 
Here we give a brief description of the components relevant to this analysis. The innermost part of CDF, a four-layer 
silicon vertex detector (SVX'), allows a precise measurement of a track's impact parameter with respect to the primary 
vertex in the plane transverse to the beam direction p]. A time projection chamber determines the position of the 
primary vertex along the beam direction. The central drift chamber, located inside a 1.4-T superconducting solenoidal 
magnet, measures the momenta of the charged particles. Outside the drift chamber there is a calorimeter, which is 
organized into electromagnetic and hadronic components, with projective towers covering the pseudo-rapidity range 
|?7| < 4.2. The muon system is located outside the calorimeter and covers the range \r]\ < 1. Events for this analysis 
were collected using a trigger which required missing transverse energy ^t > 35 GeV. ^t is the energy imbalance in 
the directions transverse to the beam direction using the raw energy deposited in calorimeter towers with \r]\ < 3.6. 

After removing events with large ^t from accelerator-induced and cosmic ray sources, we select events with two 
or three jets that have transverse energy Et > 15 GeV and \r]\ < 2 (hard jets) and no jets with 7 < Ex < 15 GeV 
and \r]\ < 3.6 (soft jets). These requirements efficiently reject ti events (which have more than 3 hard jets) and 
QCD multijet events (which have soft jets due to gluon radiation). Jets are found from calorimeter information 
using a fixed cone algorithm M with a cone radius of 0.4 in rj-(j> and jet energies are calculated using the raw energy 
deposition in calorimeter towers. The angle (j) is the angle in the plane normal to the beam direction. To reduce 
systematic effects from the trigger, we require events to have ^t > 40 GeV, and to reject events with fake missing 
energy arising from jet energy mismeasurements we require that the missing transverse energy direction is neither 
parallel to any jet(j) nor anti-parallel to the leading Et jet : A0(^t, j) > 45° and A(p{^T,ji) < 165° where the jet 
indices are ordered by decreasing Et- Moreover, to reduce the QCD background, we require the angle between the 
two leading jets to be 45° < A0(ji, J2) < 165°. We reject events with one or more identified electrons (muons) with 
Et (Pt) > 10 GeV (GeV/c). 

After applying these requirements, the data sample (which we call the pretag sample) contains 396 events. The 
largest source of background in the pretag sample is the production of M^-|-jets, where the W decays to a neutrino 
(leading to missing energy) and either an electron or muon that is not identified or a tau which decays hadronically 



The pretag sample also contains QCD multijet events where the large ^t is due to jet energy mismeasurement. 

The SVX' information is used to tag heavy-flavor jets. We associate tracks to a jet by requiring that the track 
is within a cone of 0.4 in rj-cj) space around the jet axis. We require tracks to have Pt > 1.0 GeV/c, positive 
impact parameter, and a good SVX' hit pattern. A good SVX' hit pattern consists of three or four hits in the 
SVX' detector with no hits shared by other tracks. We take the sign of a track's impact parameter to be the sign of 
the scalar product of the impact parameter and jet Et vectors. We then define the impact parameter significance 
to be the impact parameter divided by its uncertainty. For tracks originating from the primary vertex the impact 
parameter significance distribution is symmetric around zero with a shape determined by the SVX' resolution, while 
decay products of long lived objects tend to have large positive impact parameter significances. We therefore use 
the negative impact parameter significance distribution to define the detector resolution function. For each track, 
we determine the probability that the track comes from the primary vertex using this resolution function. We call 
this probability track probability. By construction, the track probability distribution is flat for tracks originating from 
the primary vertex, and peaks near zero for tracks from a secondary vertex. We combine the track probabilities for 
all tracks associated to a jet to form the jet probability {V jet) [H, the probability that all the tracks in the jet come 
from the primary vertex. The distribution of Vjet is flat for jets originating from the primary vertex by construction, 
while for bottom and charm jets it peaks near zero. Vjet is a continuous variable and the Vjet requirement is easily 
optimized for the scalar top and scalar bottom searches separately. This motivates the choice of the Vjet tagging 
algorithm over other tagging algorithms developed at CDF ||]. 

We select events for the scalar top analysis by requiring the event to have at least one taggable jet with a Vjet < 0.05. 
A taggable jet has at least two SVX' tracks as defined above. The distribution of the minimum jet probabihty ij'jll^) 
of the taggable jets in the pretag sample is shown in Figure |[ The expected background and scalar top/bottom signal 
distributions are overlaid. This requirement, chosen to optimize the expected signal significance, rejects approxiamtely 
97% of the background while its efficiency for the signal is 25%. For the scalar bottom analysis the expected signal 
significance is optimized by requiring that the event has at least one taggable jet with a P"!™ < 0.01. This requirement 
rejects approximately 99% of the background while retaining 45% of the scalar bottom signal. 

Backgrounds (other than QCD multijet events) and the expected signal are estimated using a number of Monte 
Carlo programs followed by a full CDF detector simulation. Single vector boson production and decay is simulated 
using a tree-level calculation as implemented in the VECBOS llQ] package, with HERWIG lO] routines used for 
subsequent parton hadronization. Vector boson pair production and decay is implemented in ISAJET p2[ . Top pair 
production and decay is simulated using HERWIG. Signal events are modeled using the PYTHIA flS] generator. The 
PYTHIA Monte Carlo generator includes production and decay of SUSY particles flj]. The next-to-leading order 
(NLO) cross section for the scalar quark production is calculated using the PROSPINO ||l^ program with CTEQ3M 
parton distribution functions [ p6| . Simulated events are analyzed using the same procedure as the selected data 
sample. We check the single vector boson normalization with data by reversing the lepton veto requirement in the 
pretag sample. 

We estimate the number of QCD multijet events in the tagged samples using a combination of Monte Carlo and 
data samples. We attribute the excess of data events above elcctroweak sources in the pretag data sample to QCD 
multijet sources. The total expected electroweak background in the pretag sample is 270.1 ± 76.2 which gives us an 
estimate of 125.9 ± 83.4 expected QCD multijet events in the pretag sample. We then apply a Vjet mistag matrix to 
this excess to estimate the QCD multijet background after tagging. The Vjet mistag matrix, which parameterizes the 
probability that a jet has Vjet < 0.05 as a function of jet Et and the number of SVX' tracks, is derived from data 
and verified in several control data samples. 

The systematic uncertainties on the expected number of signal events apply for both ti and bi. The NLO cross 
section for third generation scalar quarks depends weakly on other masses and parameters (~ 1%) Pq) . The dominant 
NLO uncertainties are due to the choice of QCD rcnormalization scale {fi) and the choice of parton distribution 
function. The theoretical uncertainty on the NLO scalar quark production cross section is a function of the scalar 
quark mass and ranges from 11% to 22% for the mass range 30 GeV/c^ to 150 GeV/c^. Gluon radiation from the initial 
state (ISR) or final state (FSR) partons is the largest source of systematic uncertainty. We determine its effect on our 
acceptance by turning off ISR or FSR in the signal Monte Carlo and comparing the efficiency with the default Monte 
Carlo which has ISR and FSR turned on. The combined ISR/FSR systematic uncertainty is 23%. We determine the 
jet energy systematic uncertainty, which is 10%, by varying the jet energies by ±5%. The trigger efficiency systematic 
uncertainty, which is 10%, is determined by varying the trigger efficiency curve (which is derived from data) by ±lcr 
of its fitted values. The detection efficiency estimates are derived from Monte Carlo that has exactly one primary 
vertex. The dominant effect of multiple primary vertices is to reduce the efficiency for a requirement of no extra 
jets with Et > 7 GeV and \ri\ < 3.6. We account for the loss in efficiency due to the extra jet veto by combining 
the Monte Carlo with a minimum-bias data sample (consistent with the number of primary vertices found during 



the 1994-95 Tevatron run), measuring the relative loss in efficiency of the "no extra jet" requirement, and degrading 
the signal efficiency by this factor. The efficiency scale factor due to multiple primary vertices is 0.93 ± 0.03. We 
use data samples enriched in charm(bottom) jets to determine the systematic uncertainty on the charm (bottom) 
tagging efficiency. The systematic uncertainty is 10% for both charm and bottom tagging. Including the systematic 
uncertainties due to the integrated luminosity measurement (4.1%) and finite Monte Carlo statistics (5-15%), the 
total systematic varies from 31% to 36% as a function of the squark mass. 

In the scalar top analysis we observe 11 events, which is consistent with 14.5±4.2 events expected from SM processes 
(see Table ||). We interpret the null result in the scalar top search as an excluded region in to^o-tti^^ parameter space 
using a background subtraction method [|l^. The 95% C.L. excluded region is shown in Figure |[ The maximum m^^ 
excluded is 119 GeV/c^ for m-^o = 40 GeV/c^. The maximum excluded value of the neutralino mass is 51 GeV/c^ 
which corresponds to the scalar top mass of 102 GeV/c^. The reach in rrif: is limited by the statistics, while the gap 
between the kinematic limit and the excluded region is mostly determined by the ^t cut which is effectively fixed by 
the fix trigger threshold. Also shown in Figure ^ are the results from the D0 experiment, based on 7.4 pb^^ [ [l8[ , 
and from the OPAL experiment for ^/s = 189 GeV at LEP (l|]. 

In the scalar bottom analysis five events are observed with an expected background of 5.8 ± 1.8 (see Table |). 
Similarly, we interpret the null result as an excluded region in m^o-iTif^ parameter space as shown in Figure ^. For 
m^o = 40 GeV/c^ the maximum 771^^ excluded is 146 GcV/c^ . Also plotted are the latest results from D0 pQl ^^i^d 
OPAL (19). 

In summary, we have performed a search for ti/61 in pp collisions at ^/s = 1.8 TcV using 88 pb~^ of data. We looked 
for events with significant missing energy, no high Pt lepton(s), and two or three jets. We required that at least one 
jet is consistent with originating from a heavy fiavor jet using a technique called jet probability. After applying all 
selection criteria, we observed no excess of events above Standard Model predictions, and we set 95% C.L. exclusion 
regions in the m^o-mq-^ plane. 
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Sample 


iV..p (P^T < 0.05) 


Ne.p (P™r < 0.01) 


VF*(^ e^Ue)+ > 2 jets 


0.3±0.3±0.1 


- 


VK±(-> m*/^m)+ > 2 jets 


0.9±0.5±0.3 


- 


VK±(--^ T±!/^)+ > 1 jets 


7.6±1.6±2.2 


3.0±1.0±0.9 


Z"{^ uV)+ > 2 jets 


1.2±0.4±0.4 


0.8±0.3±0.2 


tt 


0.7±0.2±0.4 


0.5±0.2±0.2 


Diboson {WW, WZ, ZZ) 


0.4±0.1±0.1 


0.2±0.1±0.1 


Total W/Z/tt/ Diboson 


11.1±1.8±3.3 


4.5±1.1±1.2 


Total QCD 


3.4±1.7 


1.3±0.7 


Total Expected from SM 


14.5±4.2 


5.8±1.8 


Total Observed 


11 


5 



TABLE I. The number of observed data and expected background events. For W/Z/tt/ Diboson, the first uncertainty is 
statistical, the second is systematic. For QCD and Total Expected from SM, the uncertainty is statistical plus systematic. 
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FIG. 1. 



The distribution of P™r 



the lowest value of Vjet for all taggable jets in an event. A requirement of 
jet ^ ^■'-"J y-j-^^! in appiicu uu Bcicuu charm (bottom) jets. Points are data, the shaded histogram is the sum of the 
predicted backgrounds, the solid line is the predicted signal for jtij =110 GeV/c^, m-a — 40 GeV/c^, and the dashed line is 



VZr < 0.05 (0.01) is applied to select charm (bottom) jets 



the predicted signal for m^ = 140 GeV/c 



40 GeV/c . The background and signal are normalized to 88 pb 
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FIG. 2. 95% CL exclusion region (shaded region) in m^o-mi plane for ti — > ex?. Also shown are results from D0 mM and 
OPAL [g^. 9{ parameterizes the mixing of the left/right scalar top gauge eigenstates to form the light/heavy mass eigenstates. 
Note that the results for both D0 and CDF are independent of d^. 
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FIG. 3. 95% CL exclusion region (shaded region) in m^o-nii plane for fei — > fexi- Also shown are the results from D0 ]2C|] 
and OPAL ||19|. 9^ parameterizes the mixing of the left/right scalar bottom gauge eigenstates to form the light/heavy mass 
eigenstates. Note that the results for both D0 and CDF are independent of 8^^. 
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